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teractions are studied. In the model, we assume that, at long distances, the two-body 
attractive potential decays as r~ a . Thermodynamic extensivity (nonextensivity) is 
observed for a > D (0 < a < D). Particular attention is payed to the liquid-gas 
critical point located, in the temperature-pressure plane, at (T C ,P C ). (T C ,P C ) are, 
in the N — > oo limit (N = number of molecules), finite for a > D and diverge for 
■ a < D (as (a — D)^ 1 for a/D — > 1 + 0). However, the variables T* = T c /N* and 

P* = P c /N* with N* = [N 1 -^ - 1]/[1 - a/D] remain finite for all a. Thus, the 
extensive and nonextensive regions become unified, as recently conjectured. These 

results should be useful for discussing gravitation and some special fluids. 
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In recent years, much attention has been paid to physical systems with microscopic long- 
range interactions. Nonextensive behavior and difficulties with the standard formulation of 
thermodynamics are since long known for this class of systems (see, for instance, ]I||§). More 



precisely, in order to have a well denned thermodynamic limit in the usual sense, a system 
must have a finite free energy per particle. In other words, an increase of N (number of 
particles) is expected to leave the free energy per particle unaltered if iV >> 1, i.e., the free 
energy Fn should be asymptotically proportional to N. However, when the effective range of 
the interactions between particles decays slowly enough with the distance, Fn increases with 
a higher power of N, hence F^/N diverges when iV — > oo. According to many text-book 
interpretations of thermodynamics, the concept itself of thermodynamic limit is ill-defined 
in such cases and the formalism becomes unphysical (see, for instance, Jl|). In some sense 
(that will become transparent later on), this is a restricted interpretation of thermodynamics 
and a possible way out from this difficulty has been recently proposed || and illustrated 
through several examples In the present effort, we shall address along these lines a 

Lennard-Jones-like fluid, with special emphasis onto its gas-liquid phase transition and the 
associated critical point. To do this we shall perform molecular dynamics for the two- and 
three-dimensional cases. 

In order to understand the thermodynamical behavior of systems with long-range inter- 
acting particles, we consider a Lennard-Jones-like model with pair interactions characterized 
by a (12, a) potential. A variety of values for a can be associated with standard interac- 
tions in models for fluids. For instance, for D = 3, a = 6 corresponds to the standard 
Lennard- Jones fluid, a = 3 essentially corresponds to dipole-dipole interactions in systems 
with dipole configurations such that the interaction is attractive, a = 2 corresponds to 
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dipole-monopole interactions (such us those relevant in the tides), a = 1 mimics the grav- 
itational and the (unscreened) Coulombic interactions, and finally a = corresponds to a 
mean field approach. For arbitrary D, the choice a = D — 2 corresponds to the isotropic 
solutions of the D-dimensional Poisson equation. Also, there are some special situations 
in which neutral, glassy, small spheres can interact through long-range potentials with low, 
and not even integer, a. Such cases might occur, through a Casimir-like effect due to the 
large fluctuations at the critical point of standard fluids, within which the little spheres are 
immersed ||. We might occasionally refer to such fluids as Casimir-like ones. For all these 
reasons, we shall focus on herein all the values of a in the interval < a < 6. 

An interesting though preliminary discussion of this class of fluids has been done by 
Grigera Our aim in the present work is to provide a more detailed description which 
will hopefully give a suitable picture of some general thermodynamic properties of systems 
with long-range (and not singular at the origin) interactions. 

The two-body potential we shall assume is given as follows: 



v(r) = 4e 
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where r is the distance between particles and a > and e > are specific Lennard- Jones-like 
parameters. For the long-tailed attractive term r~ a in potential ([!]) we consider < a < 12. 
The r -12 term describes the short range repulsive potential (due to nonbonding overlap 
between the electron clouds for the D = 3 case). The usual convention in the theory of fluids 
(for instance, standard Lennard- Jones fluids) is to define the diameter a, via v(a) = 0. Due 
to the range of the interactions, this would not seem particularly helpful here. So, we rather 
make a different convention, namely the r~ a coefficient has been chosen so that the minimal 
value w(r mm ) = e in all cases. The Lennard- Jones-like (12, a) potentials as a function of r 
are plotted in FIG [1] for typical values of a. 

A main difficulty for fluid systems with long-range (but not infinitely long) interactions is 
the absence of any exact solution. In our present case, calculations were performed with the 
molecular dynamics method for systems of N spherically-symmetric particles with periodic 
boundary conditions. Standard mean field (or van der Waals) theory for an integrable (i.e., 
a > D) potential with a cutoff r c uses corrections oc v(r)d D r. However, our main 
aim is to discuss the trend of thermodynamic quantities as a function of iV (the size of 
the system) and to provide an unified approach to the thermodynamic limit for the entire 
range (0 < a < 6). Consequently, no corrections are considered in the present treatment 
and a cutoff distance of half of the box size was applied to the interaction; we verified 
that, in the iV — > oo limit, no physical consequences seemed to emerge from the adoption 
of this computational convenience. The equations of motion were solved using the Verlet 
algorithm and the temperature was kept constant by a weak coupling of the system with 
an external thermal bath [0]. 

When starting a simulation from scratch, the initial configuration is generated by ran- 
domly distributing the molecules. Then, the molecules evolve in such a way as to achieve the 
energetically most favorable positions and velocities. Full equilibrium is assumed only when 
the pressure, the potential energy and the total energy exhibit stable values. In general, the 
computational work is relatively heavy, since for systems with long-range interactions it is 
uneasy to obtain very precise numerical data for quantities such as the critical temperature, 
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for example. Nevertheless, as we will see, the trends of the present numerical data clearly 
exhibit, in all the cases we have studied here, the conjectured thermodynamical scalings @. 

In the remainder of this paper, quantities are specified in usual reduced (dimensionless) 
form in terms of the Lennard- Jones parameters a and e. The reduced number density is po~ D 
where p = N/V; the reduced temperature is k^Tje where ks is the Boltzmann constant; 
the reduced pressure is a D P/e; finally, the reduced time is Jejm t/a, where m is the mass 
of the particles. From now on we consider a = e = m = kB = l- 

Special attention is payed to the point at which the critical isotherm (T = T c , where T c 
corresponds to the critical temperature) simultaneous has a vanishing slope (dP/dV) T = 
and an inflection point, i.e., where the curve changes from convex to concave, hence, 
(d 2 P/dV 2 ) T = in the pressure- volume plane. In fact, these properties simultaneously 
occur only at the critical point of the system. It is known that the standard Lennard- 
Jones fluid provides a handy model for testing liquid theories and for investigating such 
phenomena as melting, the liquid-vapour surface, nucleation, etc. Consistently, our interest 
in the present more general fluid stems mainly in the critical point of the liquid- vapor curve. 

The iVV^T-ensemble was used in our simulations. Naturally, strictly speaking, only at 
the thermodynamical limit we can speak of the critical temperature. Nevertheless, even for 
finite N, we can define an effective "critical" temperature T c (a; N) by using the zero slope 
and inflection point criteria just mentioned. We have obtained the values for T c (a; N) from 
simulations with systems containing iV = 27,64, 125,216 and 512 in three dimensions, and 
N = 25, 49, 100, 225 and 400 in two dimensions. 

In FIG |2|, T c (a; N) as a function of N is depicted for typical values of a for D=3 case. 
For a > D, the critical temperature shows a tendency to saturate in the limit iV —* oo, 
but the convergence becomes very slow when a approaches D from above, consistently with 
what was already discussed for the total and the potential energies of a similar fluid |§]; as 
we shall see, for a < D, the approximate critical temperature diverges with N. A suitable 
plot for the curves associated with the extensive regime (i.e., a > D) is depicted at inset in 
FIG |2|. Labels on lines correspond to several values of a at Figure (a = 1,2,3,4,5,6) and 
inset (a = 3.5,4,5,6). The D = 2 case exhibits similar trends. In fact, all our results (for 
all (a, D)) are satisfactorily fitted with the simple expression T c (a; N) » T c (a) + a a N l ~ a l D , 
where the parameters T c (a) and a a come out from the fitting. For the extensive region 
(a > D) T c (a) represents the thermodynamic limit of the critical temperature; a a < 0. 

In FIG [| we depict some values of T c (a) given by the present approach as a function of 
a. We verify that T c diverges for a < D and that, when a — > D from above, T c (a) diverges 
like l/(a — D) (see the inset in the figure). This behavior has been conjectured for generic 
systems with long-range interactions. The thermodynamic limit is not reached in the usual 
way in the nonextensive region. Nevertheless, if we appropriately scale the various relevant 
quantities, our numerical results will show that this limit is in fact as well defined as that 
corresponding to the extensive region. Let us now clarify how this happens. 

Thermodynamic quantities An like the internal energy, the free energy, Gibbs energy, 
etc., associated with systems including potentials having an attractive tail that decays as 
r~ a , do not always scale with N. Indeed, we generically have 



Consistently, we can introduce @Q] a new variable given by 

ATl-a/D _ i 

N* = — - (3) 

which, in the limit N — > oo, behaves asymptotically as l/(a/D — 1) for a/D > 1; as 
IniV for a/D = 1; and as N l ~ a ' D /(I - a/D) for < a/D < 1. For a = 0, N* ~ N, 
which precisely recovers the value frequently used in mean field approximations in order 
to renormalize the coupling constants in such a way as to make the system to (artificially) 
become extensive. It is known since several decades ( [|Kj and references therein) that 



thermodynamical extensivity imposes, in classical systems like the present one, a/D > 1. 
Our definition (3) clearly is consistent with this, but also provides the correct scaling for 
< a/D < 1, where no essential knowledge has been developed until very recently, and 
along the present lines. 

We expect the thermodynamic variables to generically include scaling with N* as recently 
conjectured ||. For instance, we expect the liniTv^oo An/(NN*) to be finite for all values 
of a, whether it is in the extensive region or in the nonextensive one. Thus, we can write 

Gn Un T Sn P Vn 

AT* AT ' \ I 



N*N N*N N* N N* N 

It can be easily verified that this type of scaling preserves the Legendre transformation 
structure of thermodynamics, even in the nonextensive region (0 < a < D) 0. These 
features are in fact included in a more general context, namely that of a recent statistical- 
mechanical formalism addressing nonextensive systems |TI|-0|. In what concerns the so 
called intensive variables (T, P, etc), we expect them to scale with iV*. For instance, the 
critical temperatures must scale as T c (a; N) ~ N*T*(a). Indeed, in FIG |(a), the critical 
temperature T c (a; N) vs iV* is exhibited for typical values of a, and a linear dependence with 
iV* is clearly observed. As before, it is possible to fit the data with the simple expression 
T c (a; N) w b a + T*(a)N*, where b a < and T*(a) are fitting parameters. Since T* is 
finite in all cases, it is to be considered as the correctly scaled expression for the critical 
point of the system. In FIG §(&), we plot T* as a function of a in both the extensive and 
nonextensive regimes. We see that T*(a) increases continuosly with a but its derivative 
dT*(a)/da possibly presents a discontinuity at a — D. The error bars illustrate the mean 
error (less than 5%) observed in the measurements of T c in our numerical simulations. 

In FIG |5](a), we plot P c (a; N)/T c (a; N) as a function of l/iV 1//3 for the D = 3 system, and 
it turns out to be once again linear. This enables a simple extrapolation (the precise finite 
size scalings associated with this problem have not been taken into account here, since they 
do not appear to add anything specially relevant for our present discussion) which provides 
P c (a)/T c (a), which of course coincides with P*(a)/T*(a), since both P and T scale with 
iV*. We notice that, as it happened with T*(a), the values are finite in both the extensive 
and nonextensive regions (i.e., Va): see FIG |5](&). The observed error (^ 10 % or less) in 
the measurement of the P c is roughly the double of that observed for T c . 

Let us finally address a more general situation, namely the functional relationships such 
as the equation of states, in the present case the funtion f(P*,p,T*) = 0, which relates 
the density p with temperature and pressure. In FIG ||(a), we illustrate the ideas with the 
D = 3 system using several sizes, namely N = 27, 64, 125, 216 at the same fixed temperature 
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T = 3.4 and a = 2. We notice, as expected, that the curves in the (P, p) plane do not 
remain finite as N increases. There is no usual thermodynamic limit and consequently no 
interesting properties can be exhibited under these conditions. The present class of models 
might produce a negative (here, unphysical) pressure P c (2; N) when N increases. Such 
behavior is commonly observed in iV-dependent isotherms in the P — p plane. However, 
and interestingly enough, a completely different picture emerges if, instead of presenting the 
results at fixed T, we do it at fixed T* = T/N*. Indeed, in FIG 0(6), we plot the curves 
corresponding to T* = 0.43 (> T*(2)) and we observe that the curves start converging as 
N increases. In fact, the convergence can be exhibited versus l/iV 1 / D (see inset), which 
allows the establishment of a well defined thermodynamic limit even for this example in the 
nonextensive region. In addition to these results, we appoint that a new family of isotherms 
is obtained in the present way, when the normalized function T*(a; N) instead of T(a; N) 
is used, but the true temperature value is always given by T(a; N) in all cases. In the 
same line, we verify that the values of P(2; N)/T(2; N) appear to converge, as desired, to 
nonnegative values when N — > oo. 

At this point, it is important to clarify that the thermodynamical variables we are an- 
alyzing are those which appear in relations such as the equation of states. The parameters 
which characterize thermodynamical equilibrium are the usual ones, i.e., the temperature T 
of the external thermal bath, its pressure P, etc. 

Summarizig, we have shown an example, a Lennard- Jones-like fluid, that can violate the 
usual scalings of thermodynamics, namely those corresponding to extensive systems, i.e., 
those associated with interactions whose range is not too long, more precisely satisfying a > 
D. We have shown how the appropriate scaling (through N* given by definition (5)) of the 
various thermodynamical quantities enables an unified picture for both short- and long-range 
interacting systems. By so doing, we have empirically established a variety of secondary 
finite-size scalings. In addition to this, we have verified that the critical temperatures, 
pressures, etc, diverge like l/(a — D) when a approaches D from above, a fact which is 
expected to be model-independent. Let us finally add that we hope to have provided enough 
arguments to show that thermodynamics and the concept of thermodynamical equilibrium 
can very well accomodate the (until now almost unexplored) nonextensive systems. The 
detailed study of more such systems (quantum, frustrated, etc) certainly is very welcome. 

One of us (S.C.) ackowledges partial support by FONDECYT, grant 3980014. The other 
one (C.T.) is thankful to Per Bak for useful remarks and warm hospitality at the Niels Bohr 
Institute where this work was concluded. Finally, CNPq and PRONEX (Brazilian Agencies) 
are also ackowledged for partial support. 
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FIGURES 



FIG. 1. Lennard-Jones-like (12, a) potentials as functions of r for a = 0, 1, 2, 4, 6. 

FIG. 2. Approximate critical temperature T c (a; N) as a function of N for the D = 3 system, 
for a = 1,2,3,4,5,6. Inset: In the extensive case (a > D), numerical data (circles) were fitted 
using T c (a; N) « T c (a) + ciaN 1 ' 01 / (solid line), for a = 3.5, 4, 5, 6. 

FIG. 3. Thermodynamic limit critical temperature T c (a) as a function of a for D = 2, 3. Inset: 
1/T c (a) versus a. For a < D the critical temperature T c diverges. 

FIG. 4. (a) The D = 3 approximate critical temperatures T c (a;N) as functions of N* for 
a = 1,2,3: Numerical data (circles) and their fitting using T c (a; N) ss b a + T c *(a)iV* (solid line); 
(b) Thermodynamic limit for the scaled critical temperature T*(a) as a function of a for the 
D = 2, 3 systems. 

FIG. 5. (a) The ratio i-" c (a; N)/T c (a; N) versus 1/iV 1 / 3 in both extensive and nonextensive 
regions for the D = 3 system; (b) The thermodynamic limit values of P c (a)/T c (a) versus a for the 
D = 2, 3 systems. 

FIG. 6. The ratio P(a)/T(a) versus the density p for a = 2 for the D = 3 system with 
iV = 27,64, 125,216 particles, at fixed T = 3.4 (a), and at fixed T* = 0.43 (b). Inset of (6): The 
ratio as a function of 1 /iV 1 / 3 for typical values of p in order to illustrate their (linear) convergence. 
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